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ABSTRACT 


A test of the equality of the parameters hy and A, of two exponential 
random variables Xy and X5 1s developed and evaluated for the case of a 
single sample from each population. The test is valid with or without 
time truncated testing although the test is less powerful under time 


truncated testing. The results are applicable to reliability studies. 
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I. INTRODUCTION 


A. PURPOSE 
The purpose of this thesis is to present a method for testing the 
equality of the scale parameters, \} and A», of two exponential random 
variables X and X» whose respective density functions are 
es > 
exp ( d1X,) r 0, xy > 0 
when time truncated testing of only one sample from each population is 
used. Time truncation consists of defining the new random variable Y 
where 
| hie he =< 
eed | 
Pee ifex = | 


B. BACKGROUND 

Lentner and Buehler [Ref.1] investigated the problem of making 
inferences about the relative magnitudes of the scale parameters of two 
untruncated exponential random variables X and Xo when sampling from 
each population was of size greater than one. In this case the sum of 
the random samples from each population is considered to be the observa- 
tion of a random variable Y; whose distribution is gamma with parameters 
A; and n-. The conditional distribution of the difference, Y, -Y5 


given the sum Y, +Y, is dependent on the parameter \y- A» . This result 


holds for any distribution from the exponential family, that is, a 
distribution of the form 


f(y) = c(Adh(y) expl7;(y) Q,(A)| 





This result was the basis for obtaining a confidence interval for the 
parameter Ay- A>: 

Lloyd and Linow [Ref. 2] present methods for obtczining confidence 
intervals for the parameter 4 in the case of a single population Y with 
an exponential distribution. These procedures utilize the chi-square 
distribution of the difference Y;,,-Y, i=], ***,n of the order statistics, 
Cay, Y, < Ys YgS +S Yn4] for observations from an exponential 
distribution. 

Cox [Ref. 3] uses the ranked observations from two exponential 
distributions and presents a modified F test for testing the hypothesis 
that both samples are from the same distribution, that is Ay= Ao. The 
power, P [the statistic falls in rejection region], of this test was 
investigated by Gehan and Thomas [Ref. 4] for the case of time truncated 
testing and small samnie sizes uStnu cuommuce: stemdatiun. Tie vest 


ee 


Statistic here is ty/t, where 


(Vql 4 sum of observed values from sample 7) 


(n--V. ) ( number of observed values from sample i} 


/ 


vee number of truncated aseeiuatrone 
The test statistic has an approxiiate F distribution with degree of 
freedom 2(n)-v) and 2(Np-Vo). Their conclusion was that Cox's modified 
F test is nearly as powerful as the standard F test. The modified F test 
is not appropriate for the case of single samples since time truncation 


yields a zero divisor. 





fee APPLICATIONS 

The results obtained herein are applicable to reliability studies of 
components assuming that the times to failure are distributed as exponen- 
tial random variables. The problem consists of choosing an alternate type 
of component to replace the type presently in use. In this case, high 
cost or large size might dictate the need to select an alternative with 
equal reliability when the items to be tested are in small supply and the 
time allocated is not sufficient for full testing. Another application 
for this procedure is the problem of making inferences when only one trial 
from each population was made, and outside considerations forced the 
termination of testing, for example, an electrical power failure during 


the initial testing of two limited production items. 





TT. DERIVATION OF PROCEDURE 


A. DISTRIBUTION OF THE SUM AND DIFFERENCE UNDER TIME TRUNCATION 
The random variable X with density function 
exp(- Ax) >0,x>0 
is transformed to the random variable Y with density function 
exp(- Ay) O<y<T 
exp(- AT) y=T 
when time truncated at time T. Using the simultaneous transformation 


Z,=¥ tat, 


eZ 
Zo=¥-Yo 
The joint density function of Z, and Lo may be found by use of Laplace 
Maeeerormations. First tre wesMey Fave of Vo. int, 2, 135 TENT Toon 


1 
WGe E exp(-At) -u(t-T) | dexp(-dt) | 


+ u(t-T) | exp(-AT)] , 
where u(t-a) is the unit step function at t=a, 
u(t-a)= 0 t<a 
= | tra 
and u (tea) is a unit mass function or Dirac delta function at t=a. 


u'(t-a) has the property 


i 


X 
fu’ (t-a)dt = 0 X4a 
~00 


to 


Now 7.2 ie So) = [E | exp (-S)Z] “S919 by definition and by 
ec 


independence of Y, and Y, this becomes 


IO 


€. (coms e—a ) 4s- +s } f, (s.-S,) 
eee ey) 1 2’ YY, 1 °2 


For this case, using standard tables of Laplace transformations, 


i a 


By multiplication of fy, (si +85) and ie 2 then inverting the 


transformation, first with respect to S, and then with respect to S53 
the result, in terms of unit step functions and unit mass functions, is 


the joint density function of Z, and Lo. This density is 


(1.1) (Agh9/2) exp | -(A,+ A )25/24(Ag-Aq )Z,/2 | O<eq<T 
IZ,1<z, 
Vea cl 
|Z51< 2T-z 
bh <2 ) N,exP LO -Ay )T-A,z, | T<z,<2T 
ee 
Gliz3)) A, exp | (Ay-Ay)T- 42, | Tez, < 27 


Z5=Z)-21 


et exp [-(, 4,7 | Zo ze 0 


when hi= a A, this becomes 
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fale). |) (¢/2) exp | Az, IS a Z|<z 
1£z1S2T, |z,|<2T~z 


iaee. 1). ) exp [ -z, | Te2z,<eT, 29=2T-2, 
ec.) ) exp [-z,d | Ts2,<2T, Z5=2)-21 
(1.4.1.) exp L-2aT ] 24=2T, 2570 


To obtain the conditional density of Z5 given Ly» the marginal density of 


Z, is obtained by integration of the joint density with respect to Zo. 


For example when T<2)<2T, 7 Cap 1S 
] 


271-2) 
(Ay Ay/2) exp[-(hy#2,}2,/24(yrA,) 2,/2) a, 


ae e | 
on \exP L(Xp- A, )T-Apz, + Ayexp | (Ay-Ap)T=Ay2, | , 


the marginal density of Z. is then 
| 


(1.5) yd, oA) [exec 4,2, )-exp(-2, 2, )f 0<z4<T 

(1.6) a -A,) | [ exp (-dp2, + XyT-AyT) 
ae + X exp | (yA TA2) 
+ Noexp| - Nz 4- (Aq) T| Tsz,<2T 

(1.7) exe +407] 1, =21 

and when \ =A \, this reduces to 

(1.5.1) (%/2)] exp (-z, )| (22, ) 06267 


(liom) (9/2) [ exp (-z,)| 2(2T-z)) 


+2 hexp(-Az,) Tsz,<2T 


Geese exp (-2AT) z,=2T. 
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Finally, the density of Z, given Z, is obtained by division of (1.1) by 
e5) ; 


(1.8) pet hy 2 xp LAZY) 29°? a er) glee 


-exp | - (Ao7 Ay )2 0<z,<T 


and dividing (1.1) by (1.6) _ 
fot) 
(1.9) (em) ( a“ exp 2 2," 2 zy 
Agex | ~2o24-Oq-2p)T | “AY exP [-hy2, 
-(No-hy )T| 


Eyal 





] 
tsz sel 


The other portions of fs 7 will not be used in the development and are 
21] 


1S 


ee ee i pelt Ye = it ndatia ; 
fmeenees. WET A.=A= A, che SeReitione! density & i) 
1 ¢ or 


’ 
a. 4 


ee) 1/22 J2]< 2)» 052467 


W/[2(2t-z,)+4/A] —— |a,\<2T-z,, T<2,< 21 
PEeurIStic ratio best, similar to a likelinood ratio test, of = Ay 


versus \)*2, using the density function f iSuautes te wiichiicaltsarer 
Vez 
rejecting the null hypothesis Hot y= Ae when R defined as 


23. @ 120) Aq A>) 


#2252122 |AytAg! 


1s less than or equal to a fixed number K and calls for accepting Hy when 


Ris greater than K. 


v4 





In the problem being considered, if such a K exists then solving 


— oe alles SO AV), < K 
(A,A,/2) exp [_(dg~Ay 2, /2- (dy 40g )2, /2] 


; , a 4 
for Z5-Zy yields a region of the form | (2, 525): IZ, zy] 9(k.2,)} 


called the critical region. Any observation falling in this region 
calls for rejecting the null hypothesis. The probability of an observa- 
tion falling in the critical region 1s called the power of the test. 

If time truncated testing had not been used and a trial 
solution of IZJ=(1-%)z, for the critical region is assumed, the power of 


the test 1s 


as — 
0 a ae, dzdz, ts a os Tle, dz5dZ, 


a\es 
yu - ~ 
7 j 


and if Ay=Ao= A, this reduces to 


2 2 
J 21XxX eXp [ -dz, | dz, ee 


When Aya the probability that an observation falls in the critical 


region is, with truncation, 


‘% (a (" pow 
£ dz dz te 7 73 dz5dz, 
0 (1-x)Z, ZZ, 2] Oey =z, We 


2T/(2-«)¢ 2T-z 2T/(2-« ) -(1-«)z, 
“| 1 dz,dz, + fi dzeaz 
J Z)lZp 2 1 f | 


T (1-x)zZ, é T 12-27 Zie> al 
pave j 
+ ies dz dz 
T Z. =eT-Z, aN. 


This 18 1] minus the value listed in Table 1. 


is 





The use of time truncation maps the region 


_ CAT all om 
A = 2912, > 2125-5 am ) Ne, Jee —, 1$21 6 a= 


2-X 


onto 


B= tee = Alea Tsz)8 | uv Gan < Iz5|< 21-2); Ts SZ, sot 


om 


Clearly B is contained in A and A is contained in the critical region, thus 
under the nuli hypothesis the probability of observing |Z] 2 (1-x)z, iis 
less than or equal ton. This is the probability of rejecting the null 
hypothesis when Ho is true. Figure 1 displays a lower bound for the 
power of this test when ~=.3 and T=21. 

If testing is not time truncated, a power function which is monotoni- 
cally increasing in the absolute value of (An-A, )2 1S available. the 


] 
test statistic is Z») given Z,, and the decision rule is: reject Hp: y=ho 


The power of the test in this case is 


: 7 
- ii dl 
a On. 4 Z| Z, “2 
exp [(A27A1) (21-21-21) /2] ~exP| = (Ap~Ay ) (2-42, 42, )/2| 
I-exp [-(Ap-Ay )z, | 


Sinh | (ho Aq) (1-00) 24/2 | 
sinh L(Ag-Ay) Z4/2 ] 


when )y=d, the power is « as may be verefied by L" Hospital's rule. 
Figure 2 displays the power function without time truncation when x=.1 and 


for selected values of Ly. 
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The probability of type II error is the probability that an observation 
falls in the compliment of the critical region, | (2, 229): Jzol <(1-x) 24 f , 


when ie is not true. This probability is 


_ ad ; a 
Z5GZ 
; -(1-x)2, i, ae 


00 zy" AT (2x) 
\ dz, dz 
2T/(2-x) 2 K2T/(2-x)-z, ve 1 


(1-4) (hy Xp)/ (2-20) hg + 25] | (z-eAgtrry | 


exp [-)y T+ do T/ (2-0) ] exp[-hoT+\y xT/(2-«) | 
(2 =x) 4 40 Ao hye ie (2e “K)AgtKA, Ay 


iO) 


+ 


Table 1 shows this probability for selected values of My odos T anda. When 


T becomes arbitrarily large, (1.10) reduces to 4 (1-0) (A 9)/ [(2-x) A, tx |- 
[(2-%)ho+ do | 


B. SUMMARY OF TEST PROCEDURE 
The procedure for testing the hypothesis Ho: Ny =Ao against the hypothesis 
Hy Ay=Ap with significance ™ 


1. Observe Yy and Y5 


2. Form 2]=¥,+Y9 and Zo=V,-Yo 

3. Compute oe 

eS Sia < (1-x)Z, accept ie 
x) 


] 
ia = (= x)Z, pegect: ti 
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II]. ESTIMATION OF Xo- Aq 


A. POINT ESTIMATION OF o> dj 


In section II, the density of Ly given Z, was derived for the case 


where neither X, nor X,are time truncated. This is the case 0<Z,<T and 


the density is given by equation 1.8. The expected value of Ly given Z, 


is then 


eee SS 


2, : 1-exph-(Xg-Ay )241 





Z | 
\ 1, (ded) exe [Og Ay) (22-2 V2] gy, 


= 7. Itexp[-(o-A; a - 2 
1-exp[- (Ap-)q) 21] 2M 


~ “1 
‘x . A 


agen re" Vo eet arts ae 
= | Ce NIC d Sarco ONIN Ce 
From this an estimate of the parameter No My may be obtained by solving 


Z5=z, coth [(do-Ay 23/2 - 2/(ho- Ay) 


fordo-)}- Numerical computations using Newton's method are appropriate 


with this equation. A typical example might be z=2; Z5=-1. 32. Then 


-1.32=2 coth (\o-Ay) = -2/ (Ag-A7) 


and the value of \o- hy which satisfies this is \o7 Ayo3: 


This result is also valid for the case T$Z)<2T and ZL 21-2, since 
neither X nor X, have been time truncated in this case. 

For the case in which X, or Xp are time truncated, a point estimate of 
Ao-Az 1S not available. But a confidence interval fordj-A7 may be obtained. 





B. CONFIDENCE INTERVALS FOR Non Ay 


A 100(1-x)% confidence interval for a parameter is an interval on the 
real line, (a,b), such that the estimator of that parameter has a 1-« 
probability of being observed in the interval (a, b) and the unknown 
parameter will lie in the interval. 

For the exponential random variable X with scale parameter A, the 
unbiased estimator for X is the observation X. A 100(1-«)% confidence 


interval for Ais 


é 


( -In(1l-«/2) , -In(«/2} | 
X 


To obtain a (1-x) 100% confidence interval for Ao-Ay when neither X nor 
X5is time truncated, first obtain separately (1-3) 100%, 4 =171-%, 
confidence intervals for \, and Aj. Then by the transformation }=Ag-Ay> 4 


confidence interval with confidence greater than or equal (1-«) 100% is 


—in¥2 4 In(i- 3/2) Sein eee ao moc 
X9 X] : Xo x) 


When either X, or ¥ are time truncated, a (1-«x) 100% confidence interval 


2 
70% 1d 1S 


_ Ink 
( 0, - 158) 


again using the transformation \=\y-A) > assuming Y;=y7<T and Vout a 


confidence interval with confidence greater than or equal (1-%) 100% is 


ee , . Ing + 1nd/2 ) 
y ip y 


where § =1-“1- 
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C. APPLICATIONS OF CONFIDENCE INTERVALS 

An ad hoc procedure for accepting the hypothesis \=A2 based on 
confidence intervals for \o-\y would be to compute the confidence interval 
based on the observed values of the random variables Y, and Yo. Then if the 
point 0 is in the interval obtained, accept H): \j=A9- This procedure 
corresponds favorably with the procedure proposed in Section II. 

To illustrate this method and its correspondence with the proposed 


procedure of Section II, consider the observations 


Y,=T=1 and Yo = .6666 
Here Z,=1.6666 and Z,=.3333. When «=.3, (1-«)Z)=1.1666 which is greater 
than Z5, hence accept Hot Apa Ag: 
For 8=1-/1-«=.165, an 83.5% confidence interval for A, is (0,-In.165)= 


(0, 1.800) and an 83.5% confidence interval for Ao is 


( (-1n.9175)/.666, (-In .0825)/.666) ) 
=a os, 3.70}! . 
Now, a confidence interval for No-Ay with confidence greater than or equal 
to 70% is 
ie leovte 3.76) 
which includes the point zero, so the hypothesis A 7A should be accepted. 


This agrees with the test result from Section II. 


2 | 





TV. APRIORI KNOWLEDGE OF A 


PeeeerHOlCE OF TEST TIME T. 

When the value of the scale parameter is known apriori, the experimentor 
may reduce the probability of type II error by a more advantageous choice 
omeenme test time l. If V5 is considered to be the standard, hence \» known 
apriori, the probability of type II error, given by equation 1.10, can be 


bounded. Specifically, the terms involving AgT may be bounded by choosing 


T appropriately. For example, if the constraint 


exp | -A,7| ey, 
is to be imposed, then T must be greater than “In ), ; 


For reliability studies, the scale parameters are usually known from 


previous data and hence can be of use for this type of test. 
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V. AREAS FOR FURTHER STUDY 


A. EXTENSION OF RESULT TO SAMPLES GREATER THAN ONE 

In order to extend the results obtained thus far to tests with more 
than one sample point from each population, it would be necessary to 
determine the joint density of the sum of the Ly; random variables and the 
sum of the Lo. random variables. Here the independence of random variables 
within populations and the independence between populations of the random 


variables makes the computations easy since 


TS : = ; is) _ H if —— (s) 


low) se a 27 ae 1 
where an equal number, m, of observations are taken from each population. 


A 7 . 7 ‘ oa! ” - ‘ AA a ¢ 
mee TANG; 2eereiding i: ceemet axtendinn the paselr wanid he ro 


> Lj. = (1x) >. 


1 1 


test 


and reject the hypothesis les S =A, if the inequality was observed. Here 
one needs the density function of this test. Without truncation, the 


density of the See Lo. given the sum mu ee is the density whose 
1 1 1 


Laplace transform is 


:| ex(- a Z rE 24 


which by independence is 


TT = ex (a4 | 7” 


1=1 


Ae 





When A=A,» this 1s the produce of m Laplace transformations of uniform 


eZ. ., Zz.) density functions, hence it is 


TT | exp(sz,.) - saat 


Vi 


j=] 2sZ.. 


i 1 / sions | 
(0) Sy 

For m larger than 2, this is difficult to invert in order to obtain the 

density function. 


Another possibility is to consider the index sets 


= {i Sas <r} 
1; ' 
De = = 
K {3 TZ) 827, 2. <21-2,,] 
an ‘ e —°/f : 
Te pa iW e Ly. a) ee) 
M = { 4 Zo i=L).-2l 


N= fi: Z,2T } 
and then consider a likelihood ratio test 


f = , 
° - ‘ . {. 
F294 9244 §224244 | AytAo) 
Using the independence of the paired observations this ratio is the product 
over the respective index sets of the densities for the region corresponding 


to the index sets. Solving this for 25. yields a critical region of the 


> aes Ev 2d, tC, d.1..\>< i, ‘ 
i€9 Sl Gex 21 I Ge] iem 27 


form 
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This suggests using the previously mentioned test 
7 


a be = (1x) > ae 
1 


but with the modification of introducing scaling factors for pairs belonging 
to index sets L and M. Such a modification might be investigated as part of 


a future study. 
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VI. CONCLUSIONS 


The test prescribed in this thesis has significance level less than 
or equal &, if time truncation is used, and significance level &, if time 
truncation is not used, when testing the equality of parameters for random 
variables with exponential distributions. As the parameter A,,- \, moves 
away from the value zero, the power function increases. 

When time truncation is not used, the power function is a monotonically 


increasing function of the absolute value of (ho-A )Z 


leas 

The test prescribed is applicable when testing with only a single 
Sample from each population and it can be used instead of Cox's modified 
F test in the case of single samples when time truncation is used. 

The use of confidence intervals described in this thesis is an 
extension of the procedure of Lentner and Buehler to the case of single 
Samples from each population. This procedure may be used in the case of 
time truncated testing; the procedure of Lentner and Buehler may not. 


However, the confidence intervals used here are not exact, while those 


of Lentner and Buehler are. 
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